Quasi-long-range order in nematics confined in random porous media 
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OS ■ We study the effect of random porous matrices on the ordering in nematic liquid crystals. The 

OS ' randomness destroys orientational long-range order and drives the liquid crystal into a glass state, 

as ■ We predict two glass phases one of which possesses quasi-long-range order. In this state the correla- 

tion length is infinite and the correlation function of the order parameter obeys a power dependence 
O ' on the distance. The small-angle light-scattering amplitude diverges but slower than in the bulk 

nematic. In the uniaxially strained porous matrices two new phases emerge. One type of strain 
■ induces an anisotropic quasi-long-range-ordered state while the other stabilizes nematic long-range 

j. ' order. 
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Quenched disorder is inevitably present even in the most pure solids. This explains a lot of phenomena, e.g. the 
residual resistance of metals. On the other hand, liquids are usually homogeneous and introducing quenched disorder 
in them requires special efforts. One of the approaches consists in pouring a liquid into a randomly interconnected 
network of pores. Such liquid-porous-matrix systems emerge in many natural and technological processes giving rise to 
a lasting scientific activity. The recent surge of interest in the field is due to a new micropore material: silica aerogel. 
Its density can be varied in a wide range up to more than 99% void volume fraction. This allows the investigation 
of both strongly and weakly confined fluids. The most interesting situation emerges in systems with many degrees of 
freedom, e.g. He-3 [0 and liquid crystals ln these substances the porous matrix not only geometrically confines 

the liquid but also induces a random orienting field that fixes the direction of the order parameter near the surface of 
the matrix. The random-field (RF) disorder is known to cause spin-glass effects J^] and such phenomena were indeed 
observed experimentally in liquid-crystal-aerogel systems. In particular, a slow glassy dynamics was reported in Refs. 
^,|). Another effect of the disorder is the suppression of the isotropic-nematic phase transition. 
In many experiments |^,^,|| the sharp isotropic-nematic transition inherent to the bulk liquid crystal was substituted 
I , by continuous ordering. The genuine phase transition was observed ^| only in highly porous aerogels. However, even 
T — I ' these systems do not have nematic long-range order (LRO) as follows from the Imry-Ma argument and the 

results of deuteron NMR measurements [R. Hence, a question arises about the nature of the new phase emerging 
below the transition. This issue was addressed in several recent numerical simulations 11 - 14 . Most of them used 
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simplified models with a discrete symmetry |11JT13| . Although the structure of the phase diagram [|1 1| |1 3|] agrees with 
the experiments the nematic LRO emerging in these models is an artifact of the discrete symmetry which allows LRO 
in three-dimensional RF systems Q. An attempt to investigate numerically a more realistic continuous model 
suggested an interesting scenario of quasi-long-range order (QLRO) that is the infinite correlation length and power 
dependences of the correlation functions on the distance. However, the results of Ref. should be taken with care 
since the Imry-Ma domain length fl5|| is comparable with the sample size used for this simulation. 

In the absence of reliable numerical results analytical approaches are especially desirable. However, the existing 
theory [ fl3[ does not extend beyond the mean field approximation. As in the other RF problems it underestimates 
£3 i fluctuations and incorrectly predicts nematic LRO. The present paper for the first time studies the low-temperature 
phase of disordered nematics beyond the mean field theory. We apply the functional renormalization group (RG) in 
• i-h , 4 — e dimensions. RG flows have different character depending on the type of ordering. LRO corresponds to the fixed 
point in which the disorder and temperature are zero. If there is only short-range order the RG flow enters a strong 
coupling regime. A non-trivial fixed point describes QLRO. It was known for a long time that QLRO is possible in 
two-dimensional pure systems with the Abelian symmetry Recently the same type of ordering was discovered 
in three-dimensional RF systems |l^,^8|. In these systems QLRO is not prohibited by the non- Abelian symmetry 
p9[ . Thus, QLRO is more common in disordered systems than clean ones. We demonstrate that QLRO emerges in 
our problem whereas LRO is absent. An experimental signature of this ordering is the divergence of the small-angle 
light-scattering cross-section (||). 

The free energy density F = Fd + F prn of the nematic in the porous matrix includes the Frank distortion energy 
p0| Fd = [Fi(divn) 2 + ^(ncurln) 2 + Ky,(n x curln) 2 ]/2, where n is the director, and the interaction F pm with the 
surface of the random matrix. The interaction tends to align the director parallel to the surface . We model the 
interaction as F pm = (hn) 2 , where h is a random vector representing the normal to the surface. This is the simplest 
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choice compatible with the equivalence of the opposite orientations of the director. Due to the universality the model 
captures all large-scale physics and allows quantitative prediction of the critical exponents in the QLRO state. The 
average amplitude of the random vector h is a measure of the disorder strength. It is a phenomenological parameter 
which depends on the pore size, anchoring energy and fractal structure of the porous matrix. The microscopic 
expression for h is discussed in Ref. Q . 

In the one-constant approximation K\ = Ki = K3 the energy F = + F pm reduces to the Hamiltonian of the 
random-anisotropy (RA) Heisenberg model |19|. Since that model possesses QLRO the same ordering is expected for 
the randomly confined nematic. However, in all ncmatics K\, K3 > K~2 and this could change the critical exponents of 
the correlation functions in the QLRO state in comparison with the random Heisenberg model. Below we demonstrate 
that this is not the case, i.e. the nematic in the porous matrix belongs to the universality class of the RA Heisenberg 
model. To get a simple idea why it occurs we first consider a two-dimensional nematic film with the director n = 



(n x , n y 



(cos fa sin fa 0) confined in the plane xy of the film in the absence of the disorder. The Frank energy is 



F d = (Ki + X 3 )(V(/)) 2 /2 + (K 3 - K 1 ){cos2fa(d x fa 2 - (d y fa 2 ]/2 + sm2fad x (j)dy<f)}. The low-temperature phase of this 
system possesses QLRO, only the term (Ki+K^CVfa 2 /2 being relevant at the large scales since (sin 20) = (cos 2fa) = 
The systematic consideration is based on the RG equations in 4 — e dimensions. Our method follows the line of Ref. 
|fl9! and is briefly described bellow. It is convenient to eliminate the disorder with the replica trick We search 
for a zero-temperature fixed point. We ascribe the scaling dimension to the director n. Then the scaling dimension 
of the temperature Ax = —2 + 0(e). Any term of the effective replica Hamiltonian containing m different replica 
indices is proportional JM to 1/T m ~ x . This allows us to show that all operators with three or more different replica 
indices are irrelevant |24j]. Hence, all relevant operators of the appropriate symmetry are included in the following 
replica Hamiltonian 
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where a, b are replica indices, a, (3 — x,y,z label the spatial coordinates, Ai = K2,X 2 = K\ — i^2,A 3 = K3 — K2, 
T is the temperature, the function R(z) describes the disorder, and the summation over the repeated indices a and 
(3 is assumed. Due to the symmetry n a «-> — n a the function R(z) is even. Below we measure the temperature in 
units of K2, and hence set Ai = 1. To define the energy in 4 — e dimensions we add to the Hamiltonian ([!]) the term 
^0 Sqfl d a n^d a n%/2 7 where a labels the coordinates in the (1 — e)-dimensional subspace, (3 = x,y,z. The stability 
conditions |2G] K\ , K3 > lead to the inequality 



A 2 , A 3 > -1. 



(2) 



We represent each replica n a (x) of the director as a combination of small-scale fields (f>f (x), i — 1,2 and a large-scale 
field n' a (x) of the unit length: 



n a (x) = n' a (x) /l - ]T(C(x)) 2 + E # ("KM, 



(3) 



where the unit vectors e"(x) are perpendicular to each other and the vector n' a (x). The fields fa change at small 
scales a < r < L, where a is the molecule size, L S> a. The field n' changes at the scales r > L. The RG procedure 
consists in integrating out the small-scale fields fa and the rescaling such that the effective Hamiltonian of the field 
n' would have the structure (Q) with new constants. The rescaling is defined in such a way that Ai = 1 remains 
unchanged. The RG equations in the first order in e = 4 — D read 



dT 
d\nL ~ 

where the constant 



{D - 2)T + (1 - A 3 )C T; ^ - -A 2 (l + A 3 )Q,; ■ ^ 



dR(z = l)/dz 



dhxL 



-(3A 3 + A 3 — \2)C < f ) , 



87rVAo(l + A 3 ) 
describes the fluctuations of the small-scale fields 



1 + A 2 



(4) 



(5) 



h\) = {fa 2 ) = C^HL/a). 



(6) 
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We omit the RG equations for Ao and R(z) since their structure is irrelevant below. Eqs. have the only fixed 
point compatible with the stability conditions In this fixed point T = A 2 = A3 = and Eq. (|l]) reduces to the 
Hamiltonian of the RA Heisenberg model which thus describes the large-distance physics of the randomly confined 
nematic. Since that model possesses QLRO in its low-temperature phase for weak disorder, QLRO is also possible in 
confined nematics. For strong disorder or high temperature the ordering disappears. Thus, there are three phases: 
the high-temperature isotropic phase and two low-temperature glass phases with and without QLRO. In both glass 
phases the local orientation of the director is fixed by the random potential. As discussed below the disorder driven 
transition between the glass phases is related with topological defects. 

The large-scale correlations of the director lead to strong small-angle light scattering. We determine its intensity 
in the limit of the weak optical anisotropy, i.e. assuming that in the dielectric tensor e Q( 3 = e±S a p + t a n a nfj the 
anisotropic term e a C ei. In this case the scattering cross-section can be found with the Born approximation. The 
scattering cross-section with the change of the wave vector by q is given by the expression |2(J 

(r(q) = \u 2 /(4irc 2 )i a e a p(q)f b \ 2 , (7) 

where u) is the light frequency, i and f are the unit vectors specifying the initial and final polarizations, e Q ^(q) is the 
Fourier transform of the dielectric tensor. Hence, cr(q) ~ (Qa l s(q)Qce l 9( — <j))i where Q a p = n a np — 6 al 3/3 is the order 
parameter and the angular brackets denote the disorder and thermal averages. In contrast to the bulk nematic the 
scattering is caused not by the thermal fluctuations but by the frozen configuration of the director. The cross-section 
<x(q) is proportional to the Fourier transform of the correlation function G(r) = (Q a p(0)Q a p(r)) . In the QLRO state 
this correlator obeys a power dependence on the distance G(r) ~ r~ n . To calculate the exponent rj we decompose 
Q a p into small-scale and large-scale parts with Eq. (|J) and average over the small-scale fluctuations with Eq. (||) : 

(Q a0 (O)Q<*f3(r)U = K(tH(0)(l - + ^ e 2 Q (0)4(0)(<^-} - 6 a p/3} x 

i ij 

K(r)4(r)(l - +5>Ur)4(r)<^) -6 aB /3} = Q^(0)Q' Q/3 (r)[l - GC^nL/a], (8) 

i ij 

where Q' a g = n' a n 'p ~ S a p/3, (...) denotes the average over the fluctuations of tf>, and the relation (<fii<fij) ~ which is 
valid in the RA Heisenberg fixed point is used. The constant C<f, = 0.309e Eq. (^) is the same as in the fixed point of 
the RA Heisenberg model [[I9| . The exponent rj can be found with the iterative use of Eq. (||) at each RG step until 
the scale L = r is reached. At the scale r the values of the renormalized director field n' are the same at the points 
and r. Hence, Q' a g(0)Q' a g(r) ~ 1 and r~ v ~ (1 — QC^lnL/a) , where K = ln(r/a)/ln(L/o) is the number of the 
RG steps. Thus, 77 = 6(7^. The small-angle scattering cross-section is given by the expression 

<r(q) ~ q- D+r > = <T 4+2 - 9e . (9) 

The uniaxial stress modifies the large-distance behavior. The compression along the z-axis can be described by 
adding to the Hamiltonian the term Fg = An 2 , where A > 0, since the deformation tends to make the pore surfaces 
parallel to the xy plane and hence favors the planar configuration of the director. The uniaxial stretch is described 
by Fs = An 2 z with a negative A. In both cases A is proportional to the deformation. The effect of the electric field 
is analogous to the effect of the stress but the sign of the electric energy |^0| F e = — e a (nE) 2 /87r is fixed for a given 
substance. The RG flow is unstable with respect to the perturbation F s and new regimes emerge at the scale R = R c 
at which the renormalized A(R) ~ 1. The critical length R c can be found analogously to the correlation length of the 
RA Heisenberg model in the uniform magnetic field Q. At small A the result is R c — \A\- 1/( - 2 ~ 2C ^ = |yl|-0-5-o.i5e_ 
The stretched system is long-range-ordered at the scales R > R c . The nematic order parameter can be calculated 
analogously to the magnetization of the RA Heisenberg model in the uniform magnetic field and is given by the 
formula Q = {n a n fJ - 5 aB /Z) - Rc * ~ |A| 46e . LRO can also be achieved by applying an arbitrarily weak external 
magnetic field to the confined nematic since the magnetic contribution to the energy |2(J F m — — Xa(nH) 2 /2 has 
the same structure as the energy related with the uniaxial stretch. A more interesting situation emerges under the 
compression. The director averaged over a scale R > R c is confined in the xy-plane. The system is thus described by 
the RA XY model. It possesses QLRO but the critical exponents are different from the exponents of the Heisenberg 
model. Thus, at the scale R c the cross-over from one QLRO state to another occurs. Using the RA XY fixed 



point found in Ref. |19| and repeating the derivation of Eq. (g) one finds the Born light-scattering cross-section for 
q < 1/ R c : cr(q) ~ q- 4 + £ ( 1 + Tr / 9 ). i n the RA XY regime the cross-section (^) is anisotropic: the small-angle scattering 
is suppressed, if the incident or scattered light is polarized along the compression direction. 
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Our RG procedure is based on the decomposition (|3j) which makes sense only if the director change is slow at the 
microscopic scale a. This condition is broken in the core of a topological defect. However, the topological defects are 
irrelevant at small e = 4 — D for weak disorder. This can be understood from the consideration of the contribution of 
the disclination loops and the pairs of point defects of size I 3> a to the RG equations at the scale I. After averaging 
over the small-scale fluctuations the size I of the topological excitations plays the role of the ultra-violet cut-off. The 
renormalized temperature is small: T(l) <C 1. Hence, the thermal fluctuations are irrelevant. The disorder-induced 
term i?(n a rifc) ~ e in the renormalized replica Hamiltonian ([!]) is of order (ft 4 ), where the random vector h describes 
the (renormalized) random anisotropy F pm = (hn) 2 , (...) denotes the average over the realizations of the disorder. 
The excitation energy can be compensated by the interaction with the disorder only in the positions where h ~ 1. 
The concentration of such positions is exponentially small ~ exp(— 1/e). Hence, the defects produce the corrections 
of order exp(— 1/e) to the RG equations and do not modify the results of the paper qualitatively. The concentration 
of the topological excitations of size I is of order l~ D exp(— 1/e). The above discussion is valid, if the disorder is 
weak. In the case of the strong disorder the topological defects are present at the microscopic scale a and QLRO is 
absent. Thus, the topological defects drive the system into another glass state in which the orientation of the director 
is determined only by the local random potential. 

Recently it was suggested that the RF model (Q) describes nematic elastomers |25|| , However, it is unclear if this 
model and hence the above results are applicable to nematic elastomers because of elasticity-mediated non-local 
interactions in these substances p6|. The non-local dipole interactions are present also in the amorphous magnets 
which could be described by the RA Heisenberg model in the absence of the dipole forces. The effect of the long-range 
interactions on the stability of QLRO is an open question. Such forces are absent in randomly confined nematics which 
provide a genuine realization of the RA Heisenberg model and can be used for an experimental test of the possibility 
of QLRO in the non-Abelian systems. In conclusion, we have demonstrated that weakly disordered nematics possess 
QLRO in their low-temperature phase. The uniaxial stretch and the external magnetic field stabilize LRO while the 
uniaxial compression drives the liquid crystal into another QLRO state. 

The author thanks A. Kamenev, V.V. Lebedev, V.P. Mineev and V. Steinberg for useful discussions and R. Whitney 
for the critical reading of the manuscript. 
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